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Abstract 

The ATLAS detector at the Large Hadron Collider at CERN is used to search for the decay of a scalar boson to a pair of long-lived 
particles, neutral under the Standard Model gauge group, in 20.3 fb 1 of data collected in proton-proton collisions at y/s = 8 TeV. 
This search is sensitive to long-lived particles that decay to Standard Model particles producing jets at the outer edge of the ATLAS 
electromagnetic calorimeter or inside the hadronic calorimeter. No significant excess of events is observed. Limits are reported on 
the product of the scalar boson production cross section times branching ratio into long-lived neutral particles as a function of the 
proper lifetime of the particles. Limits are reported for boson masses from 100 GeV to 900 GeV, and a long-lived neutral particle 
mass from 10 GeV to 150 GeV. 
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1. Introduction 

The discovery of the Higgs boson OHD by the ATLAS and 
CMS experiments 13 ED in 2012 identified the last piece of 
the highly successful Standard Model (SM). Subsequent mea¬ 
surements of the Higgs boson branching ratios and couplings, 
while consistent with the SM expectations, allow for a substan¬ 
tial branching ratio to exotic particles. This letter describes a 
search for decays of the Higgs boson and other scalar bosons to 
non-SM states that in turn decay to SM particles. 

A number of extensions of the SM involve a hidden sector 
that is weakly coupled to the SM, with the two connected via 
a communicator particle. This letter considers models contain¬ 
ing a hidden sector with a confining gauge interaction that is 
otherwise invisible to the SM. The communicator is chosen 
to be a SM-sector scalar boson, <f> jbjfT The communicator 
mixes with a hidden-sector scalar boson, <5>h s , which decays 
into detectable SM particles. This search considers commu¬ 
nicator masses between 100 GeV and 900 GeV. A 4> mass close 
to the mass of the discovered Higgs boson is included to search 
for exotic decays of the Higgs boson. 

A Hidden Valley (HV) model (8JO is used as the benchmark 
model. The lightest HV particles form an isospin triplet of pseu¬ 
doscalar particles which are called valley pions (n v ) because 
of their similarity to the SM triplet. The 7T V are pair-produced 
(Ohs ~> 7Tv Jlv) and each decays to a pair of SM fermions. The 
7T V possess Yukawa couplings to fermions and therefore pref¬ 
erentially decay to accessible heavy fermions, primarily bb, cc 
and T + T~. 

The lifetime of the 7T V is unconstrained and could be quite 
long. A long-lived /r v can result in signatures that traditional 
searches fail to detect. If a 7t v decays in the inner detector or 
muon spectrometer, it can be reconstructed as a displaced ver¬ 
tex. However, standard vertex-finding algorithms oia are not 
likely to reconstruct it without modification. Likewise, a 7 r v de- 
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cay deep inside the calorimeter is reconstructed as a jet with 
an unusual energy signature that most traditional searches re¬ 
ject as having poor data quality. This search focusses on fi¬ 
nal states where both 7T V decay in the hadronic calorimeter or 
near the outer edge of the electromagnetic calorimeter. Each 
heavy fermion pair from a 7 r v decay is reconstructed as a single 
calorimeter jet with three characteristic properties: a narrow ra¬ 
dius, no tracks from charged particles matched to the jet, and 
little or no energy deposited in the electromagnetic calorimeter. 

Scalar boson masses ranging from 100 GeV to 900 GeV 
are considered in addition to the Higgs boson’s mass (gener¬ 
ated at mu = 126 GeV) and 7 r v masses between 10 GeV and 
150 GeV are studied. Other searches for pairs of displaced ver¬ 
tices generated by pair-produced neutral, long-lived particles 
were performed in ATLAS ITTI and CMS fl2l at the LHC and 
in DO lfl3l and CDF lfl4l at the Tevatron. The Tevatron exper¬ 
iments and CMS searched for displaced vertices in their track¬ 
ing system only, which results in a corresponding proper decay 
length range of a few meters. CMS also looked at the multi¬ 
lepton decay channel, another possible decay of HV particles. 
The previous ATLAS analysis, based on 7 TeV data, used the 
muon spectrometer and is sensitive to proper decay lengths be¬ 
tween 0.5 m and 27 m, depending on the benchmark model. No 
evidence of physics beyond the SM was found. 


2. The ATLAS detector 

The ATLAS detector [0! is a multi-purpose detector at the 
LHC, consisting of several sub-detectors. From the interaction 
point (IP) outwards there are an inner detector (ID), electro¬ 
magnetic and hadronic calorimeters, and a muon spectrometer 
(MS). The ID, immersed in a 2 T axial magnetic field, provides 
tracking and vertex information for charged particles within the 
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pseudorapidity Q(rj) region |rj| < 2.5. It consists of three dif¬ 
ferent tracking detectors. From small radii outwards, these are 
a silicon pixel detector, a silicon microstrip tracker (SCT) and a 
transition radiation tracker (TRT). 

The calorimeter provides coverage over the range |rj | < 4.9. 
It consists of a lead/liquid-argon electromagnetic calorimeter 
(ECal) at smaller radii surrounded by a hadronic calorimeter 
(HCal) at larger radii comprising a steel and scintillator-tile 
system in the barrel region (| rj | < 1.7) and a liquid-argon sys¬ 
tem with copper absorbers in the endcaps (1.5 < hi < 3.2). 
The ECal spans the range 1.5 m < r < 2.0 m in the barrel 
and 3.6 m < M < 4.25 m in the endcaps. The HCal covers 
2.25 m < r < 4.25 m in the barrel and 4.3 m < |z| < 6.05 m 
in the endcaps. There is also a forward calorimeter (FCal), 
with coverage between 3.1 < hi < 4.9, which uses copper ab¬ 
sorbers in the first layer, and tungsten absorbers in the second 
and third layers, and liquid-argon as the active medium in all 
layers. Muon identification and momentum measurement are 
provided by the MS, which extends to |tjI = 2.7. It consists 
of a three-layer system of gas-filled precision-tracking cham¬ 
bers. The region |tj| < 2.4 is also covered by separate trigger 
chambers. 

A sequential three-level trigger system selects events to be 
recorded for offline analysis. The first level consists of custom 
hardware that implements selection on jets, electrons, photons, 
T leptons, muons, and missing transverse momentum or large 
total transverse energy. The second and third levels add charged 
particle track finding and refine the first-level selections with 
progressively more detailed algorithms. 

3. Data and simulation samples 

All data used in this analysis were collected during the 2012 
LHC proton-proton run at a centre-of-mass energy of 8 TeV. 
After data quality requirements are applied, the sample corre¬ 
sponds to an integrated luminosity of 20.3 fb~ 1 . The HV Monte 
Carlo (MC) samples are generated with PYTHIA 8.165 Ifl6l 
and the PDF MSTW2008 fl7l to simulate gluon fusion gg —> <t> 
production and the Ohs decay Oj ls —» K V K V for different <t> and 
7T V masses (Table [T}. O masses below 300 GeV are consid¬ 
ered low-mass samples and the rest are considered high-mass 
samples. The n v lifetime is fixed in each sample to ensure de¬ 
cays throughout the ATLAS detector. The <t> is simulated in 
PYTHIA by replacing the Higgs boson with the O and having 
the <t> decay to 7r v 100% of the time. The O samples are pro¬ 
duced with cross sections calculated at next-to-next-to-leading- 
logarithmic accuracy in QCD processes and at next-to-leading- 
order in electro-weak processes assuming the O at each mass 
has the same properties as the SM Higgs boson |fl8l . After 
generation the events are passed through a detailed simulation 


1 ATLAS uses a right-handed coordinate system with its origin at the nom- 
inal interaction point (IP) in the centre of the detector and the z-axis along the 
beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the 
y-axis points upward. Cylindrical coordinates (r, (j)) are used in the transverse 
plane, 0 being the azimuthal angle around the beam pipe. The pseudorapidity 
is defined in terms of the polar angle 6 as 7] = — In tan(0/2). 


of the detector response with GEANT4 lfl9l [20 1 and the same 
reconstruction algorithms as are used on the data. GEANT4 
needed no modification to simulate the signal as all decay parti¬ 
cles are SM particles. All MC samples are reweighted to repro¬ 
duce the number of interactions per bunch crossing observed in 
the data. 

Table 1: The <3> mass or Higgs boson mass, <3> gluon fusion production cross 
section, and 7T V mass of each benchmark Hidden Valley model generated. The 
cross-sections are based on the assumption in the benchmarked model that the 
<t> boson production mechanism is the same as the Higgs boson production 
mechanism. The decay branching ratios of the 7i v as a function of the 7T V mass 
are listed in the second table as determined in the simulation samples. 


m H 

a 

7T V Mass 

[GeV] 

[pb] 

[GeV ] 

126 

19.0 

10, 25, 40 


4> Mass 
[GeV] 

a 

[pb] 

7T V Mass 
[GeV] 

100 

29.7 

10, 25 

140 

15.4 

10, 20, 40 

300 

3.59 

50 

600 

0.52 

50, 150 

900 

0.06 

50, 150 


i r v Mass 
[GeV] 

BR bb 
[%] 

BR T + T 
[%] 

BR cc 
[%] 

10 

70.0 

16.4 

13.4 

20 

86.3 

8.0 

5.6 

25 

86.6 

8.1 

5.3 

40 

86.5 

8.5 

5.0 

50 

86.2 

8.8 

4.9 

150 

84.8 

10.2 

4.8 


4. Trigger and event selection 

Candidate events are collected using a dedicated trigger, 
called the CalRatio trigger BTI . which looks specifically for 
long-lived neutral particles that decay near the outer radius of 
the ECal or within the HCal. The trigger is tuned to look for 
events containing at least one narrow jet with little energy de¬ 
posited in the ECal and no charged tracks pointing towards 
the jet. At the first level the trigger selects only narrow jets 
by requiring at least 40 GeV of transverse energy (Ej) in the 
calorimeter in a 0.2 x 0.2 (A rj x A0) region using topological 
jets mm, in contrast to the default algorithm in which the 
energy in a 0.4 x 0.4 region is summed. The 40 GeV Ej thresh¬ 
old requirement is fully efficient at an offline jet E T of 60 GeV. 
To select jets with a high fraction of their energy in the HCal 
the second level of the trigger requires these narrow jets to have 
logjQ^H/EEM) > 1-2, where Eh/Eem is the ratio of the energy 
deposited in the HCal (Eh) to the energy deposited in the ECal 
(Eem)- The trigger also requires no tracks with p\- > 1 GeV in 
the region 0.2 x 0.2 (A rj x A0) around the jet axis. The third 
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level of the trigger uses the slower but more accurate anti-A:, al¬ 
gorithm lf23j with R = 0.4 to reconstruct the jet and requires the 
jet to have a minimum of 35 GeV of transverse energy. 

The probability (e^) for a single 7T V to fire the trigger in sim¬ 
ulated events is shown in Fig. |T| for the (a) barrel and (b) end- 
cap region of the calorimeter in several different signal samples. 
The average probability for the low (high) scalar boson masses 
is about 20% (55%) for 7T V decays occurring at radii between 
2.0 m and 3.5 m in the barrel, and about 6% (30%) for 7T V de¬ 
cays with |z| between 4.0 m and 5.5 m in the endcaps. The 
turn-on takes place before the inner edge of the HCAL as the 
log 10 (£ , H/£EM) cut allows for a small amount of energy in the 
ECal. The probability decreases towards the outer region of the 
HCal where too much of the energy escapes the HCal to pass 
the jet Ej requirement. The efficiency is lower in the endcaps 
because events tend to not satisfy the isolation criteria due to 
the increased occupancy from extra collision events in the same 
bunch crossing as a hard-scatter interaction (pile-up). 
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Events also contain a reconstructed primary vertex with at 
least three tracks with pj > 1 GeV. Events are rejected if any 
reconstructed jets show evidence of being caused by a beam- 
halo interaction ETTl . A missing transverse momentum require¬ 
ment, £™ ss < 50 GeV, is applied to reject non-collision events, 
such as cosmic rays or beam-halo interactions. 

In the offline selection, jets are reconstructed with an anti-A:, 
algorithm with R = 0.4, starting from calorimeter energy clus¬ 
ters calibrated using the local cluster weighting method l24l . 
Jets are then calibrated using an energy- and r;-dependent 
simulation-based calibration scheme. Jets are rejected if they 
do not satisfy the standard ATLAS good-jet criteria with the 
exception of requirements that reject jets with small electro¬ 
magnetic energy fraction (EMF) [12511 . At least one jet must 
have fired the CalRatio trigger. The jet matching the trigger 
must pass an Ej > 60 GeV requirement while a second jet must 
satisfy an Ej > 40 GeV requirement. If more than one jet fired 
the CalRatio trigger then only the leading jet is required to have 
Ej > 60 GeV. 

Individually, all jets must satisfy |r)| < 2.5, have 
log id(Lji/£|-; m) > 1.2, and have no good tracks in the ID with 
pj > 1 GeV in a region A R < 0.20 centred on the jet axis. A 
good track must have at least two hits in the pixel detector 
and a total of at least nine hits in the pixel and SCT detectors. 
Fig.[2ja) compares the distribution of the number of good tracks 
associated with each jet in the multi-jet sample (described in the 
next section) with that in jets resulting from simulated 7T V de¬ 
cays in the HCAL or ID. Fig.[2jb) makes the same comparison 
for the distribution of log 10 (£ , H/£EM) of each jet. The multi-jet 
data was gathered using a prescaled, single-jet trigger with a 15 
GeV requirement. 

Jets caused by cosmic rays and beam-halo interactions are 
often out-of-time. The jet timing is calculated by making an 
energy-weighted average of the timing for each cell in the jet. 
Each cell is defined to have a time of 0 ns if its energy is 
recorded at a time consistent with the arrival of a J3 = 1 par¬ 
ticle from the IP. The timing of each jet is required to satisfy 
— 1 < t < 5 ns. This cut will impact the efficiency for low p 
7 r v . Due to the requirement of a high-Ex jet in this analysis the 
P -distribution is peaked near 1 for low mass <t> samples. For 
the high mass <t> samples the difference between iihj and 
results in a large boost for the 7 r v at the generated lifetimes. As 
a result, the inefficiency introduced by the timing cut is at worst 
1.5% for the considered samples. 

The analysis requires that exactly two jets satisfy these re¬ 
quirements. The second jet requirement significantly reduces 
the SM multi-jet background contribution. Table [3] lists the fi¬ 
nal number of expected events in each signal MC sample. The 
final number of events selected in data is 24. 

5. Background estimation 

The largest contribution to the expected background comes 
from SM multi-jet events. Cosmic-ray interactions contribute 


Fig. 1: The probability (£ Ky ) for a single 7T V to pass the trigger as a function of 
the 7i y (a) radial decay length in the barrel and (b) the z position of the decay 
vertex in the endcaps for several and 7T V masses. 
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Fig. 2: Distribution of (a) the number of good tracks (bracks) with pj > 1 GeV 
and A R ■v 0.2 around the jet axis and (b) the distiibution of jet lo^ jq 
with jet 1 77 1 < 2.5, f> \ > 40 GeV. The dashed histogram is for n., jets decaying 
in the hadronic calorimeter, and the dotted histogram is for 7T V jets decaying in 
the ID. Both are from the in// = 126 GeV, m Iv = 10 GeV sample. The filled 
histogram is the multi-jet data sample used to evaluate the multi-jet contribution 
to the background. Events are required to satisfy Ej' ss < 50 GeV. 


at a much lower level, and beam-halo interactions make a neg¬ 
ligible contribution. 

To estimate the multi-jet background contribution, a multi¬ 
jet data sample is used to derive the probability that a jet passes 
the trigger and analysis selection. To obtain a raw background 
prediction, these jet probabilities are applied to a data sample 
that represents the multi-jet background before application of 
jet-level analysis selection. A correction to account for two-jet 
correlations is applied to this raw prediction to yield the final 
multi-jet background estimate. 

The multi-jet data sample contains events that pass single¬ 
jet triggers with an Ej threshold of 15 GeV or higher. These 
triggers were prescaled in 2012 and their effective luminosi¬ 
ties range between 14.8 nb^'and 454.1 nb _1 . The dataset is 
representative of the full 20.3 tb 1 of data collected in 2012 
and contains events from throughout the data collection pe¬ 
riod. The events are required to pass the analysis E™ lss require¬ 
ment and have at least two back-to-back (A (j) > 2.0) jets with 


E\ > 40 GeV and —2.5 < 77 < 2.5. One jet is required to sat¬ 
isfy the modified ATLAS good-jet criteria used by the analysis. 
The second is used to measure two probabilities: one, called P, 
for a jet to pass the trigger and the Ej > 60 GeV jet require¬ 
ment and the other, called Q, for a jet to satisfy the require¬ 
ment Ei > 40 GeV. For both P and Q the jet must also pass 
the log 10 (£ , H/£EM)j track isolation, and all other analysis jet 
selection requirements including the modified ATLAS good-jet 
criteria. The probabilities are determined as a function of jet 
Ej and 77. The jet E \ and 77 dependence is calculated indepen¬ 
dently because the sample is not large. A systematic uncertainty 
to account for any potential correlation is included in the anal¬ 
ysis. To calculate this systematic uncertainty the change in the 
mean E\ as a function of 77 and the change in the mean 77 as 
a function of Et is measured in the multi-jet data sample. The 
maximum variation is 2% for both E-\ and 77 . The Ej or rj of 
each jet are systematically shifted by this amount as P and Q 
are recalculated. The new P and Q distributions are used to 
estimate the multi-jet background as described below, and the 
maximum variation in the result (6%) is used as the systematic 
uncertainty. 

Binned fits of the probabilities as a function of Ej are made 
to a Landau function and an exponential function for P and Q, 
respectively. The Ej requirement is ignored when fitting to al¬ 
low the curve to best match the distribution’s shape. The fit 
errors are propagated through to the systematic in the multi-jet 
background. The 77 dependence is strongly correlated with the 
distribution of material in the calorimeters and cannot be well 
described by any simple functional form. Thus the probability 
is obtained directly from the distribution. The P(Ej) parameter- 
isation is additionally split into leading jet and sub-leading jet 
samples because the probability is different for the two types of 
jets. This effect is also present for Q; however, it is accounted 
for by the correction for jet correlations discussed below. Plots 
of P(Ej) and Q(E^) are shown in Fig. [ 3 ] The peak present 
in P(Ej) is the result of the trigger turn-on for the full trigger 
chain. The trigger jets are dominated by leading jets, and so 
dominated by the leading jet P. 

The probability P is verified using the CalRatio-triggered 
data. The CalRatio-triggered events are required to pass the 
same event selection used to derive the single-jet probabilities 
as well as the requirements for calculating P. The CalRatio- 
triggered data contains 501 387 events that fired the unprescaled 
CalRatio trigger and passed the required selection, and the 
single-jet probabilities predict 513 000 ±94000 (statistical er¬ 
ror only) events. 

To calculate the raw multi-jet background prediction the 
probabilities P and Q are applied to jets in events selected by 
the 15 GeV single-jet trigger. These single-jet probabilities are 
combined into an event probability using a combinatoric cal¬ 
culation that requires at least one jet in the event to fire the 
trigger and exactly two jets to pass all the jet selection crite¬ 
ria. The event probability is scaled to account for single-jet 
trigger prescales, yielding a weight for each event. The sum of 
all weights in the data sample yields a raw background predic¬ 
tion of 13.2 ± 2.9 (statistical ® systematic error) events. The 
uncertainty is dominated by the small number of jets firing the 
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Fig. 3: The probability, P, that a jet from the multi-jet data sample passes 
the trigger and all jet requirements including the Ej > 60 GeV requirement is 
shown in (a). A Landau function is fitted to the leading and sub-leading jet 
distributions separately (solid and dashed lines). The probability, Q, to pass all 
jet requirements including the Ej > 40 GeV requirement as a function of jet Ej 
is shown in (b). An exponential function is fitted to the distribution (solid line). 
The Ej requirement is ignored when fitting to allow the curve to best match the 
shape of the data, but is used explicitly when P and Q are applied to a jet. The 
multi-jet data was gathered using a range of prescaled, single-jet triggers. 

CalRatio trigger in multi-jet events. 

In multi-jet events the log ]q (£h /-Eem ) and track isolation 
values of one jet are correlated with those of the second jet. 
If an event contains one jet of high log^En/EEM), the second 
jet is more likely to have high log m (/f|-i/£’i:Mj as well. Like¬ 
wise, if one jet has no tracks associated with it, the other is 
more likely to have no associated tracks as well. The single-jet 
probabilities above ignore this correlation because each is cal¬ 
culated independently of the log^En/EEM) and ;; tra ck of other 
jets in the event. As a result, Q is lower than if it were cal¬ 
culated only in events with an accompanying low «track- high 
1 °gio( £, H/£ , EM) jet. 

A scale factor to account for the correlation is calculated 
from the multi-jet data sample and the CalRatio-triggered data 
sample by examining numbers of events in regions in the 


log [q (£h/Lem ) and number-of-tracks («track) plane that are out¬ 
side the signal region (log l0 (£H/£EM) >1-2 and n t rack = 0). 
The logjojEn/EEM) binning is chosen such that binning is uni¬ 
form in EMFrl A range from 0 to 7 was used for u t rack- The 
regions outside the signal region are expected to have very little 
signal contamination. 

In each region the ratio of the number of events observed 
in the CalRatio-triggered data to the raw prediction is calcu¬ 
lated. Two series of ratios are calculated, one as a function of 
\og w (E H /E EM ) and one as a function of n lrac k. To determine 
the trend in the ratio as a function of log 10 (£’H/£’EM) the n t rack 
requirement is held constant: a jet is required to have 5 or 6 
tracks. The ratio is then determined for several non-overlapping 
ranges of log^En/EEM)- The same procedure is used for /) trac k 
by requiring jets to have 0.55 < EMF < 0.65. Because the ratio 
is taken with respect to the observed data, this ratio will correct 
for any normalisation errors in P or Q. 

Both sets of ratios are fitted to allow extrapolation into the 
signal region. The product of the two ratios in the signal region 
yields a scale factor to correct for the correlation between jets. 
A systematic error is added to account for the assumption that 
the two ratios are uncorrelated. The calculated scale factor is 
1.8 ± 0.5. The uncertainty on the scale factor is due to the 
limited sample size. 

To verify the procedure eight other bins on the 
logjo^n/FEM) and track plane were chosen and the full 
background prediction method was applied. Because signal 
contamination is negligible outside of the signal region, the 
predicted number of events can be directly compared to the 
number of events in the same log 10 (£’H/FEM)-«track region 
in the CalRatio-triggered data. In all cases the prediction is 
consistent with data to within one standard deviation. 

The final multi-jet prediction is 23.2 ± 8.0 (statistical © sys¬ 
tematic error) events in the signal region. The uncertainty is 
dominated by the statistical uncertainty, which is in turn domi¬ 
nated by the small number of jets matching the CalRatio trigger 
in the multi-jet data sample. The systematic contribution comes 
from the correlation between Ej and r) as well as from the in¬ 
clusion of a requirement on A (j) (not used in the signal selection) 
in the determination of P and Q. 

Particles from a cosmic-ray shower may pass through and 
deposit energy in the calorimeter without passing through the 
ID. These energy deposits can be reconstructed as trackless jets. 
The overall contribution to the expected background is reduced 
by the jet-timing and £™ ss requirements. 

The cosmic-ray background was studied using a trigger sim¬ 
ilar to the CalRatio trigger, but active only during an empty 
crossing. Each proton beam is divided up into buckets, most of 
which are filled with protons. An empty crossing occurs when 
an empty bucket in each beam coincides in the centre of the 
detector, and five buckets on either side in each beam are also 
empty. Data gathered from these empty crossings are used to 
study backgrounds that are not beam related. 

The analysis selection, with the exception of the jet-timing 
requirement and the good-vertex requirement, are applied to all 


3 log 10 (£ H /£EM) = log 10 ((l — EMF)/EMF). 
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events triggered in empty crossings. The —1 < t < 5 ns timing 
requirement is removed to retain more events to give a more 
accurate determination of the background. A simple scaling can 
be used to predict the expected cosmic-ray event rate within the 
timing window because the arrival time of cosmic-ray muons 
is uniformly distributed. It is found that about 5% of cosmic- 
ray events firing the trigger and containing two jets are events 
where both jets satisfy the —1 < t < 5 ns requirement. 

Two additional corrections are applied to determine the fi¬ 
nal background prediction due to cosmic-ray events. The first 
accounts for the different live-times of the triggers. The num¬ 
ber of empty crossings is 2.9 times smaller than the number 
used to collect the full data of 20.3 tb 1 . The second correc¬ 
tion weights each event to account for soft tracks due to pile-up 
and underlying-event effects that would have caused the jet to 
fail the track isolation requirement had it occurred in a collision 
environment. To determine the weights a trigger that selects 
random collision events is used to determine the probability as 
a function of T] that a track with pj > 1 GeV is present in a 
A R <0.2 cone anywhere in the detector as a function of t]. 
This probability is applied to each jet in each event to deter¬ 
mine an event weight. The event weights range from 0.55 to 
0.63. Combining all the corrections results in a predicted num¬ 
ber of cosmic-ray events of 0.3 ±0.2 (statistical error). 

Another possible background contribution comes from a 
beam-halo muon that undergoes bremsstrahlung in the HCal. 
Two selection criteria reduce this type of background. A jet¬ 
timing requirement is imposed because most of the jets pro¬ 
duced by beam-halo interactions are not coincident in time with 
jets from pp interactions. In addition, events are rejected when 
track segments in the endcap muon chambers, from the entering 
beam-halo muon, align in (j) with a jet. These two requirements 
reduce the background considerably with no discernible effect 
on the signal. 

Unpaired isolated crossings, i.e. crossings where only pro¬ 
tons from a single beam are present and at least three buckets 
on either side of the empty beam’s bucket are also empty, can be 
used to study beam-halo events. To estimate this background, 
artificial events are created by sampling two jets from a col¬ 
lection of jets passing both a CalRatio trigger active only dur¬ 
ing unpaired isolated crossings and the leading jet requirements 
from unpaired isolated crossings. All possible pairs of jets are 
used and the /7 l | rHSS < 50 GeV requirement is applied to each 
constructed event. The number of jets passing the jet analysis 
selection and the fraction of constructed events satisfying the 
Ej' ss requirement are combined to estimate the background. 
This method, which also accounts for cosmic-ray muon con¬ 
tamination, predicts 0.07 ± 0.07 events. The large uncertainty 
is due primarily to the small number of jets passing all required 
cuts. 

Backgrounds from combinations of these non-beam interac¬ 
tions, i.e. a beam-halo jet plus a multi-jet, or a beam-halo jet 
plus a jet due to a cosmic-ray muon, were found to be negligi¬ 
ble. 


6. Systematic uncertainties 

Table [2] presents a summary of systematic uncertainties as¬ 
sociated with the signal sample. The overall uncertainty, taken 
as the sum in quadrature of all positive and negative contribu¬ 
tions respectively, is listed in the last column. The MC signal 
samples’ statistical uncertainty is shown in Table [3] and it is 
accounted for in the statistical analysis. The overall normalisa¬ 
tion uncertainty of the integrated luminosity is 2.8% obtained 
following the same methodology as that detailed in Ref. lf26l 
from a preliminary calibration of the luminosity scale derived 
from beam-separation scans performed in November 2012. The 
uncertainties on the Higgs boson production cross-sections at 
y/s = 8 TeV, which are equal to the uncertainties on the <t> pro¬ 
duction cross-sections, are about 10% ED- 

The uncertainty on the signal MC samples due to parton 
distribution functions (PDF) is calculated by reweighting each 
event using three different PDF sets (MSTW2008nlo68cl 11771 , 
CT10 f27l . and NNPDF2.3 ll28l ) and their associated error sets. 
The RMS change in acceptance for the error sets of each PDF 
is calculated and combined with the difference in acceptances 
for each of the three PDFs. 

Pile-up primarily affects the acceptance by adding extra 
tracks and degrading the track isolation of a jet. All MC sam¬ 
ples are reweighted to reproduce the observed distribution of 
the number of interactions per bunch crossing in the data. To 
determine if pile-up is simulated properly in the MC samples, 
a direct comparison of data and MC multi-jet samples is per¬ 
formed. The jet Ej , EMF, t], 0, associated tracks and timing 
distributions as a function of the mean number of pile-up in¬ 
teractions are compared in data and MC simulation. A 10% 
systematic uncertainty is assigned to the acceptance covers all 
the observed differences. 

The jet energy scale (JES) uncertainty is evaluated as a func¬ 
tion of the jet EMF and rj, following the same strategy used in 
the in situ jet energy intercalibration ll24l . The JES is rederived 
for low EMF jets for this analysis. The relative jet calorimeter 
response is studied by balancing the transverse momenta of di¬ 
jets. The systematic uncertainty is obtained by comparing the 
/?'[ -balance in data to the /?|-balance in MC samples. This dif¬ 
ference is used to calculate a difference in the JES in data and 
MC simulation and is propagated to the signal MC samples to 
get a systematic uncertainty on the acceptance. This study also 
provides a useful performance comparison between data and 
MC jets that resemble signal jets. 

The Ej' ss uncertainty accounts for variations in missing 
transverse momentum scale and resolution l29l . The timing 
systematic accounts for mismodeling of jet timing between 
MC and data. Both of these uncertainties were determined by 
smearing the associated cut to determine the impact on the ac¬ 
ceptance. 

The simulation of the trigger is verified by comparing the 
performance of the trigger in data with the performance in MC 
simulation on the same multi-jet sample used to evaluate the 
multi-jet background. Each trigger requirement is studied indi¬ 
vidually: the jet Ej, the log | {| (±7 1 /A- iiM) 5 and the track isolation. 
Each requirement is adjusted to make the performance match in 
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data and MC events, and the resulting differences in acceptance 
from the nominal acceptance, for each requirement, are added 
in quadrature to determine the systematic uncertainty. 

The simulation of initial state radiation (ISR) cannot be di¬ 
rectly verified because it is difficult to uniquely identify ISR 
jets in data (30). An incorrect ISR rate in the simulation im¬ 
pacts the acceptance by altering the number of jets in the event 
and by altering the boost of the <f> boson. Each of these is stud¬ 
ied independently. The ISR jet population is altered event by 
event so that the number of ISR jets is halved or doubled (jets 
in MC samples are labelled as containing ISR if they contain a 
gluon with p-\ > 2 GeV). The population of 7T V jets is not altered 
by this process, but an added ISR jet may overlap with one of 
the 7T V jets. The effect of a boost caused by an ISR jet is studied 
by exploiting the correlation between the 7i v jet Ej and the <t> 
boost. From Ref. (30l . the <t> pj spectrum has an uncertainty of 
5%, which directly correlates with a 5% uncertainty in the 7T V 
jet energy. To calculate the systematic uncertainty associated 
with the boost, the p \ of ISR jets is conservatively varied by 
5% and the change in acceptance is observed. 

The changes in acceptance from both sources of ISR uncer¬ 
tainty are taken as correlated systematic errors and added to get 
the total systematic for ISR simulation. 

An incorrect simulation of final state radiation (FSR) has a 
negligible effect on the analysis’ acceptance. FSR can occur 
in a prompt or displaced jet. But even if displaced, the extra 
jet cannot degrade track isolation or deposit extra energy in the 
ECal if the 7i v has decayed in the HCal. 


7. Results and exclusion limits 

The global acceptance of the selected event topology in the 
signal MC samples is a function of in, j>, m n .. and the proper de¬ 
cay length of the 71 v . At a proper decay length of 1.5 m the 
acceptance ranges from 0.07% to 0.61%. The main efficiency 
loss is due to the low probability that both 7i v decay inside the 
calorimeter. High mass samples suffer further efficiency loss 
due to the E'™ ss requirement. Table [ 3 ] lists the expected num¬ 
ber of events from all signal MC samples and the background 
expectation in 20.3 fb -1 . The mu — 126 GeV mass samples 
use the SM Higgs boson cross-sections of Osm = 19.0 pb for 
the gluon fusion process: other production modes are ignored. 
The number of events observed in data, 24, is also shown for 
comparison. No excess of events is observed since the expected 
background is 23.5 ±8.0. The CL S method EH is used to derive 
an upper limit on the <7(<1>) xBR(<l> —> 7T V 7i v ). A profile likeli¬ 
hood ratio is used as the test statistic and a frequentest calculator 
is used to generate toy data. The likelihood includes a Poisson 
probability term describing the total number of observed events. 
Systematic uncertainties are incorporated as nuisance parame¬ 
ters through their effect on the mean of the Poisson functions 
and through convolution with their assumed Gaussian distribu¬ 
tions. The number of expected events in signal MC samples, to¬ 
gether with the estimate of expected background, the observed 
collision events and all the systematic uncertainties are pro¬ 
vided as input for computing the CL S value, which represents 


Table 2: Summary of systematic uncertainties for the <t> and Higgs boson pro¬ 
duction cross-section, jet energy scale, trigger, missing transverse momentum, 
and the requirement on jet timing as a percentage of the signal yield. Systematic 
errors that have common values across samples are not listed (pile-up at 10%, 
ISR at j^i' 2 %, and PDF at 2.1%). The last column reports the total systematic 
uncertainty (including the luminosity and common systematic errors). 


Sample 

He 

JES 

Trigger 

miss 

Time 

Total 

m H , m, t v 

[%] 

[%] 

[%] 

[%] 

Cut 

[%] 

[GeV] 





[%] 


126, 10 

+ 10.4 

- 10.4 

+ 2.2 

-2.1 

±1.1 

+ 5.5 

- 2.4 

+ 1.6 

-6.6 

+ 16.4 
-16.7 

126, 25 

+ 10.4 

- 10.4 

+1.5 

-1.6 

±1.3 

+ 3.1 

- 1.8 

+0.8 

-3.3 

+ 15.6 
-15.5 

126, 40 

+ 10.4 

- 10.4 

+2.6 

-6.2 

±1.1 

+7.7 

-4.6 

+1.9 

-5.9 

+ 18.2 
-16.9 


Sample 

T> a 

JES 

Trigger 

miss 

Time 

Total 


m 

[%] 

[%] 

[%] 

Cut 

[%] 

[GeV] 





[%] 


100, 10 

+11.1 

-10.6 

+ 2.3 

- 4.0 

±0.1 

+ 4.6 

- 3.4 

+ 2.7 

-9.5 

+ 16.7 
-18.5 

100, 25 

+ 11.1 

- 10.6 

+ 5.5 

- 3.7 

±1.2 

+ 3.4 

- 2.5 

+ 1.7 
-0.7 

+ 17.0 
-15.8 

140, 10 

+ 10.1 

- 10.3 

+ 0.6 

- 1.1 

±0.5 

+ 4.0 

- 5.6 

+ 1.9 
-6.6 

+ 15.6 
-17.2 

140, 20 

+ 10.1 

- 10.3 

+ 1.2 

- 1.6 

±1.0 

+ 4.0 

- 3.9 

+0.4 

-5.0 

+ 15.5 
-16.2 

140, 40 

+ 10.1 

- 10.3 

+ 1.3 

- 1.6 

±1.5 

+ 6.3 

- 4.6 

+ 1.8 
-2.4 

+ 16.5 
-15.8 

300, 50 

+ 9.6 

- 10.0 

+ 0.1 

- 0.3 

±0.3 

+ 9.0 

-1.4 

+0.5 

-3.0 

+ 13.9 
-13.3 

600, 50 

+ 11.2 

- 10.1 

+ 0.0 

- 0.1 

±0.2 

+11.7 

-11.3 

+2.2 

-4.4 

+ 17.0 
-16.2 

600, 150 

+ 11.2 

- 10.1 

+ 0.2 

- 0.2 

±0.3 

+11.5 

-10.2 

+2.7 

-5.3 

+ 17.5 
-15.1 

900, 50 

+ 12.8 

- 11.5 

+ 0.0 

- 0.1 

±0.1 

+12.6 

-9.7 

+1.0 

-3.7 

+ 18.5 
-15.9 

900, 150 

+ 12.8 

- 11.5 

+ 0.2 

- 0.3 

±0.2 

+11.8 

-10.9 

+0.9 

-2.5 

+ 18.1 
-16.3 


the probability for the given observation to be compatible with 
the signal + background hypothesis. 

The acceptance is a function of the <f> mass, the 7r v mass and 
7T V proper decay length. To extrapolate to the number of ex¬ 
pected events at different proper decay lengths, a large sample 
of 7t v decays is generated in a range from 0 to 50 m and an ef¬ 
ficiency map as a function of 7 r v boost is used to determine the 
efficiency at each decay length. The resulting efficiencies are 
then converted into the final number of expected events shown 
in Fig. [4] Finally, Fig.[5]and Fig.[6]show the observed limit dis¬ 
tribution for the three 126 GeV Higgs samples and for the other 
<t> samples respectively. The derived 95% confidence level (CL) 
excluded ranges of proper decay length are listed in Table[4]for 
the mu = 126 GeV samples, under the alternative assumptions 
of a 30% BR or a 10% BR for H -P 7l v 7t v . 
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Table 3: Summary of expected number of signal events, expected background 
present in the data sample, and the observed number of events in 20.3 fb -1 . 
The global acceptance is also given. The error on the signal samples is statis¬ 
tical only, the error on the expected background is statistical ® systematic. All 
results are normalised for a proper decay length of the 7T V of 1.5 m. A 100% 
branching ratio for <3>h s —> 7r v 7r v is assumed. 



(a) 



(b) 

Fig. 4: Number of events expected to pass the analysis selection in 20.3 fb -1 as 
a function of the 7T V proper decay length for (a) three low-mass datasets and (b) 
three high-mass datasets. 100% branching ratios for <J>h s —> 7r v 7r y are assumed. 


Sample 

Expected 

Global 

{m H , m k, [GeV ]) 

yields 

acceptance (%) 

126, 10 

536 ± 23 

0.139=1=0.006 

126, 25 

941 ± 44 

0.244=1=0.011 

126, 40 

365 = 1 = 31 

0.095 = 1 = 0.008 


Sample 

Expected 

Global 

(m H , m [GeV ]) 

yields 

acceptance (%) 

100, 10 

440 ± 29 

0.073 ± 0.005 

100, 25 

424 = 1 = 37 

0.070 ± 0.006 

140, 10 

525 = 1 = 20 

0.168 ±0.006 

140, 20 

900 = 1 = 37 

0.287 ± 0.012 

140, 40 

641 = 1 = 30 

0.205 ±0.010 

300, 50 

444=1= 11 

0.609 ± 0.015 

600, 50 

35=1= 1 

0.330 ±0.010 

600, 150 

41 = 1 = 2 

0.386 ±0.015 

900, 50 

3.5 ±0.1 

0.304 ±0.011 

900, 150 

4.6 ± 0.2 

0.397 ± 0.016 


Background 

Expected events 

SM Multi-jets 

23.2 ± 8.0 

Cosmic rays 

0.3 ± 0.2 

Total Expected Background 

23.5 ±8.0 

Data 

24 



Fig. 5: Observed 95% CL limits on ct/<7sm for /w// — 126 GeV as a function of 
the 7T V proper decay length: the solid line is for m Kv = 10 GeV, the short-dashed 
line is for m ny/ = 25 GeV, the long-dashed line is for m Kw = 40 GeV. The <7sm 
is taken to be 19.0 pb. The horizontal solid line corresponds to BR = 30% and 
the horizontal dashed line to BR = 10%. 
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(a) Observed 95% CL limits for m<$ = 100 and 140 GeV. 



(b) Observed 95% CL limits for m<j» = 300, 600, and 900 GeV. 

Fig. 6: Observed 95% CL limits on ax BR [pb] for (a) low and (b) high-mass 
<f> samples. 


Table 4: Ranges of 7T V proper decay lengths excluded at 95% CL assuming a 
30% and a 10% BR for a mu = 126 GeV. 


MC sample 

Excluded range 

Excluded range 

'«//. m, r v 

30% BR H -4 7t v 7i v 

10% BR H -A n v 7t v 

[GeV] 

[m] 

[m] 

126, 10 

0.10-6.08 

0.14-3.13 

126, 25 

0.30-14.99 

0.41-7.57 

126, 40 

0.68-18.50 

1.03-8.32 


8. Summary and conclusions 

A search for the decay of a scalar boson in the mass range 
from 100 GeV to 900 GeV, including a search for an exotic 
decay of the Higgs boson, to a pair of long-lived neutral par¬ 
ticles decaying in the ATLAS hadronic calorimeter has been 
presented. The analysis is based on 20.3 fb _1 of pp collisions 
at y/s = 8 TeV collected in 2012 by the ATLAS experiment at 
the LHC. 


No significant excess of events is observed over the back¬ 
ground estimate. Limits are set on the 7T V proper decay lengths 
for different scalar boson and 7 r v mass combinations. For a SM 
Higgs decaying to 7i v proper decay lengths between 0.10 m and 
18.50 m assuming a 30% BR are ruled out, and between 0.14 
m and 8.32 m assuming a BR of 10%. Results for low mass <J> 
(100 GeV and 140 GeV) and high mass (300 GeV, 600 GeV, 
and 900 GeV) have also been presented as a function of proper 
decay length. 
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